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Abstract 

We provide a probabilistic analysis of the d-dimcnsional banker algorithm when transition 
probabilities may depend on time and space. The transition probabilities evolve, as time goes 
by, along the trajectory of an ergodic Markovian environment, whereas the spatial parameter 
just acts on long runs. Our model complements the one considered by Guillotin-Plantard and 
Schott where transitions are governed by a dynamical system, and appears as a new (small) 
step towards more general time and space dependent protocols. 

Our analysis relies on well-known results in stochastic homogenization theory and investigates 
the asymptotic behaviour of the rescaled algorithm as the total amount of resource available 
for allocation tends to the infinity. In the two dimensional setting, we manage to exhibit three 
different possible regimes for the deadlock time of the limit system. To the best of our knowledge, 
the way we distinguish these regimes is completely new. 

Keywords: distributed algorithms; random environment; stochastic homogenization; reflected 
diffusion. 
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1 Introduction 

Many real-world phenomena involve time and space dependency. Think about option pricing: the 
behavior of traders is not the same when stock markets are opening and a few minutes before 
closure. Multi-agents problems, dam management problems are typical examples where (often 
random) decisions have to be made under time and space constraints. 

In computer science such problems are usually called resource sharing problems. Their statement 
is as follows (see Maier [20], Maier and Schott |21| ): 
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Consider the interaction of q independent processes Pi, ... , P q , each with its own memory needs. 
The processes are allowed to allocate and deallocate r different, non-substitutable resources (types 
of memory): Ri, . . . , R r . Resource limitations and resource exhaustions are defined as follows. 
At any time s, process Pi is assumed to have allocated some quantity yj (s) of resource Rj (both 
time and resource usage are taken to be discrete, so that s € N and yf (s) G N). Process Pi is 
assumed to have some maximum need rriij of resource Rj, so that 

< yj(s) < (1.1) 

for all s. The numbers m« may be infinite; if finite, it is a hard limit that the process Pi never 
attempts to exceed. The resources Rj are limited, so that 

X>?(*)<A, (1.2) 
i=i 

for Aj — 1 the total amount of resource Rj available for allocation. Resource exhaustion occurs 
when some process Pi issues an unfulfillable request for a quantity of some resource Rj. Here 
'unfulfillable' means that fulfilling the request would violate one of the inequalities 1)1. 2JI . 

The state space Q of the memory allocation system is the subset of N gr determined by (jl.l)) 
and (|1.2|) . This polyhedral state space is familiar: it is used in the banker algorithm for deadlock 
avoidance. Most treatments of deadlocks (see Habermann JJ) assume that processes request and 
release resources in a mechanical way: a process Pi requests increasing amounts of each resource Rj 
until the corresponding goal rriij is reached, then releases resource units until y\ = 0, and repeats 
(the r different goals of the process need not be reached simultaneously, of course). This is a 
powerful assumption: it facilitates a classification of system states into 'safe' and 'unsafe' states, 
the latter being those which can lead to deadlock. However it is an idealization. Assume that 
regardless of the system state, each process Pi with < y\ < rriij can issue either an allocation or 
deallocation request for resource Rj . The probabilities of the different sorts of request may depend 
on the current state vector (yj). In other words the state of the storage allocation system is taken 
as a function of time to be a finite-state Markov chain; this is an alternative approach which goes 
back at least as far as Ellis [5j. 

The goal is to estimate the amount of time r until memory exhaustion occurs, if initially the r types 
of resources are completely unallocated: y\ = for all i, j. The consequences of expanding the 
resource limits Aj — 1 and the per-process maximum needs rriij (if finite) on the expected time to 
exhaustion are particularly interesting for practical applications. There has been little work on the 
exhaustion of shared memory, or on 'multidimensional' exhaustion, where one of a number of in- 
equivalent resources becomes exhausted. D. E. Knuth [15] : A. Yao [32]) P- Flajolet [B], G. Louchard 
and R. Schott |18j . have provided combinatorial or probabilistic analysis of some resource sharing 
problems under the assumption that transition probabilities are constant. Maier provided a large 
deviation analysis of colliding stacks for the more difficult case in which the transition probabilities 
are nontrivially state-dependent. More recently N. Guillotin-Plantard and R. Schott t 10. analysed 
a model of exhaustion of shared resources where allocation and deallocation requests are modelled 
by time dependent dynamic random variables. 

In this paper, we analyse such problems when the probability transitions are time and space de- 
pendent. We incorporate in the transitions of our model, the influence of an environment which 



2 



evolves randomly albeit in a stationary manner on the long run. Our analysis relies on well-known 
results in stochastic homogenization. We focus on the regime when memory exhaustion is caused 
by normal fluctuations (as in many of the above references), but not by large deviations (as in Maier 
20 ) . This paper can be viewed as a (small) step towards the analysis of protocols where decisions 
are time and space dependent random variables. In addition to incorporating space dependence, 
we improve on the results of Guillotin-Plantard and Schott JO] in covering larger perturbations of 
time- homogeneous random walks. 

The organization is as follows: Section 21 presents an example of resource sharing algorithm usually 
called banker algorithm. Precise statements of the main results are given in Section G3 Proofs are 
postponed to the three following parts of the paper: we investigate the asymptotic behaviour of 
the algorithm in Sections El and El and we discuss in Section El several different regimes for the limit 
process. As a conclusion, we indicate further conceivable extensions in Section [7| 

2 Example: the Banker Algorithm 

2.1 Practical Description of the Model in Dimension Two 

We consider a simple distributed algorithm which involves only two customers C\ and C2 sharing 
a fixed quantity of a given resource M (money). There are fixed upper bounds mi and m,2 on how 
much of the resource each of the customers is allowed to use at any time. The banker decides to give 
to the customer Cj, i = 1,2, the required resource units only if the remaining units are sufficient 
in order to fulfill the requirements of Cj, j = 1,2; j ^ i. This situation is modelled (see Figure 1) 
by a random walk in a rectangle with a broken corner , i.e., 



where the last constraint generates the broken corner. The random walk is reflected on the sides 
parallel to the axes and is absorbed on the sloping side. 



{(xi, X2) € Z , < x\ < mi, < X2 < m2, x\ + X2 < A} 




Figure 1: Banker algorithm 
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Again the hitting place and the hitting time of the absorbing boundary are the parameters of 
interest. This algorithm has been analysed by Louchard et al.|17 | ITH1 119j in the model where the 
transitions are constant and by Guillotin-Plantard and Schott ^U] when the transitions are some 
time dependent dynamic random variables. 

Our goal here is to analyse such distributed algorithms in an ergodic markovian environment where 
the transition probabilities are time and space dependent, in a sense that we specify below. The 
environment takes care of purely unpredictable behaviors in addition to periodic and quasi-periodic 
ones. 

For the sake of simplicity, we assume from now on that m\ = mi = m and that A writes A = Am,, 
with A G [1,2). The size of the problem is then given by the parameter m, which therefore plays 
a crucial role in the time scale. For large m's, we will prove approximations for the deadlock time 
by taking advantage of averaging properties of the environment. 

2.2 Probabilistic Modelization 

We generalize the situation described in Subsection 12.11 and assume that there are d customers 
instead of two. We then modelize the resources of these customers by a Markov chain with values 
in the broken corner {(xi, . . . , x^) G ^ d '■ < X{ < m, 1 < i < d, x\ + • • • + x d < Am}, m G N* and 
A G [l,d). 

Walk without Boundary Conditions. 

Define first the transitions of the chain without taking care of the boundary conditions. Recall 
to this end that the transition matrix of a time-space homogeneous random walk to the nearest 
neighbours in Z rf , d > 1, reduces to a probability p(-) on the 2d directions of the discrete grid, 

V = {ei, -e x , ...,e d , -e d } , 

where (ej)i<j<d denotes the canonical basis of R d . For given i G {1, . . . , d} and u G V, p(u) simply 
denotes the probability of going from the current position x to x + u. 

Assume for a moment that the transition matrices are space homogeneous but depend on time 
through some environment (evolving with time) with values in a finite space E, N = \E\. We 
then need to consider, not a single transition probability, but a family •), . . . , p(iV, •) of N 
probabilities on V. If the environment at time n is in the state i G E, then the transition of the 
Markov chain at that time is governed by the probability p(i, •). 

When the environment is given by a stochastic process (£ n )n>o on E, the jumps (J n )n>o of the 
random walk are such that, for every u G V: 

F{j n+1 =u\^ n J ] = P (Cn,u), 
with T*k J = cr{£ , ...,£ n ,Jo,---, Jn}- 

From now on, we assume that the environment (£ n )n>i is a time-homogeneous Markov chain on 
E, and we denote by P its transition matrix, P(k,£) = P(£ n +i = £\£ n = k) for k,£ G E. Then, the 
couple (£ n , J n ) n >i is itself a time-homogeneous Markov chain on the product space SxZ 6 * governed 
by the following transition: 

\/k G E, Vu G V, P{e„+i = k,J n+1 = u\^ J } = P(Ca,k)p(Cn,u). 
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Define the position of the walker in Z d : 



S = 0, Vn> 0, S n+1 = S n + J n+1 . (2.1) 

In view of the applications mentioned above, the model is not fine enough. For this reason, we also 
assume that the steps (J n )n>i depend on the walker position in the following way: for all k E E 
and u £ V, 

P{e„+i = k, J n+1 = u\^ J ] = P(£ n , k)p(C n , S n /m, u), (2.2) 

where m denotes a large integer that refers to the size of the box in Subsection 12.11 and, for each 
k £ E and y £ R d , p(k,y,-) a probability on V. Note that the random walk (S n = Sn)n>0 
depends on the parameter m. Nevertheless, for simplicity we will often forget the dependence on 
m in our notations. In other words, (£ n , S n ) n >o defines a Markov chain with rates: 

Vk EE, Vu £ V, P{£n+i = k, S n+1 = u + S n \^' S } = P(t n ,k)p(£ n , S n /m,u), 

where jFl' 5 = a{£ , . . . , £ n , S , . . . , S n }. 
Walk with Reflection Conditions. 

Our original problem with reflection on the hyperplanes = 0, i £ {l,...,d} and yi = m, i £ 
{l,...,d} follows from a slight correction of the former one. The underlying reflected walk 
(Rn)n>o (also denoted by (-Rn"^)n>o to specify the dependence on m) satisfies, with = 

a {£,0, ■ ■ ■ ,£n, -RO) • • • ) Rn}'- 

Vk £E,Vu£ V, P{C„+i = k,R n+1 =u + R n \^ R ) = P^ n ,k)q(^ n ,R n /m,u), 

where q denotes the kernel: 

Vk £ E, Vy £ [0, l] d , VI £ {1, . . . , d}, q(k, y, ±e e ) = p{k, y, ±e t ) if < y t < 1. 

If y t = 1, then q(k,y,e e ) = and q(k,y,-e e ) = p(k,y,e e ) + p(k,y,-e e ). If y e = 0, then 
q(k, y, -en) = and q(k, y, e e ) = p(k, y, e e ) + p(k, y, -eg). 

The deadlock time of the banker algorithm is then given by = inf{n > 0, Yle=i(R>n e f) — 
This also writes: T< m ) = inf{n > 0, R n £ mF }, with F = {y £ [0, 2] d , J2i<d \Vi ~ 1| < - A)} 
(see Figure 13 below). This latter form will be useful in the proof of the main results. 

2.3 Main Assumptions 

In formula (|2.2[) . the division by m indicates that the dependence of the transition kernel on the 
position of the walker takes place at scale m. For large m, the space dependence is mild, since we 
will assume all through the paper the following smoothness property: 

Assumption (A.l). The function p is twice continuously differentiable with respect to y with 
bounded derivatives. In particular, there exists a constant K > such that: 

Vx£E,u£V, (y,y')£(R d ) 2 , \p(x,y,u)-p(x,y',u)\<K\y-y'\. 

It is then readily seen that the transition kernel 1)2. 2|l weakly depends on the space position of the 
walker: a step of the walker modifies only slightly the transition kernel. 
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Figure 2: Absorption set. 

We also assume the environment to be ergodic and to fulfill the so-called central limit theorem for 
Markov chains. We thus impose the following sufficient conditions: 

Assumption (A. 2). The matrix P is irreducible on E. It is then well known that P admits a 
unique invariant probability measure, denoted by fj,. 

Assumption (A. 3) The matrix P is aperiodic. In particular, it satisfies the Doeblin condition: 

3m > 0, 3r]> 0, V(M) G £ 2 > (-P m )(M) > V- 

Assumption (A. 4) Define for k G E andy G R d , g(k,y) = ^2 ue y[p{k, y, u)u\ (i.e., g{k,y) matches 
the expectation of the measure p(k, y, •) ) and assume that g(-, y), seen as a function from E to M. d , 
is centered with respect to the measure fj,. 

Let us briefly describe the role of each of these conditions: 

1. Thanks to Assumption (A. 2), the Markov chain (£ n )n>o satisfies the ergodic theorem for 
Markov processes. 

2. The general central limit theorem for Markov chains with finite state space follows from the 
Doeblin condition, given in Assumption (A. 3). 

3. Assumption (A. 4) permits to apply the previous central limit theorem to the function g. 
2.4 Frequently Used Notations 

For a square integrable martingale M, (M) denotes the bracket of M (do not mix up with the 
Euclidean scalar product, which is denoted by (x,y) for x,y G R d ). 

We also denote by B(R + , R d ) the path space of right continuous and left limited functions from R + 
into R d . A function in D(R + ,M d ) is then said to be "cad-lag" for the French acronym "continue a 
droite-limite a gauche" . 



6 



3 Main Results 



3.1 Asymptotic Behaviour of the Walk without Reflection 

In light of Assumptions (A. 1-4), we expect the global effect of the environment process (£ n )n>o to 
reduce for large time to a deterministic one. To this end, we view the process (S^ = m~ l SJ^2 t \ )t>o 
clS cl random element in the space B(R+,R ): 

Theorem 3.1 The process (St)t>o converges in law in B(R + ,R rf ) endowed with the Skorokhod 
topology towards the (unique) solution of the martingale problem starting from the origin at time 
zero and associated to the operator: 

1 d d 2 d - d 

i,j=l 1 i=l yl 

The limit coefficients a and b are given by: 

VyGR d , a(y) = [ [a + gv* + vg* - 2gg t ](i,y)dn(i), 
Je 

Kv)= I [(y y v -V y g)g](i,y)dfi(i), 
Je 

where V stands for the gradient and a(i, y) denotes the second order moment matrix of the measure 

V« G E, Vy G R d , a{i,y) = ^ [Wp(i, y, u)] , 

uev 

and v(i,y) denotes the well-defined sum: 

Vi eE,\/ y e R d , v(i, y) = J2Y, i pn ^ fiM> 2/)] ■ 

n>0 j£E 

Moreover, for every y G R d , the lowest eigenvalue of a(y) is greater or equal than the lowest 
eigenvalue of the matrix: 

/ [a{i,y) - gg l {i,y)]dii{i). 
Je 

In particular, if the covariance matrices of the measures (p(i, y, '))ieE,yeR d are uniformly elliptic, 
the matrices (a(y)) y ^d are also uniformly elliptic. 

The existence of v, V y v and V y g will be detailed in the sequel of the paper. 

We can denote by cx(x), for x G R rf , the nonnegative symmetric square root of a{x). Then, the SDE 
(stochastic differential equation) associated to C writes: 

dX t = b{X t )dt + <j{X t )dB t , (3.1) 

where {Bt)t>0 denotes a (f-dimensional Brownian motion on a given filtered probability space 
(fi, (^t)t>o,P). Under Assumptions (A. 1-4), (|3.1|) is uniquely strongly solvable. 

Of course, if p(k,y, •) does not depend on y (as in the beginning of Subsection 12.2(1 . then the drift 
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b vanishes in Theorem 13 . 1 1 since 6 = and the matrix- valued function a is constant so that X is a 
non-standard Brownian motion. 

Theorem 13.11 appears actually as an homogenization property: on the long run, the time inho- 
mogeneous rescaled walk behaves like an homogeneous diffusion. Underlying techniques to 
establish the asymptotic behaviour of (S^ )t>o are well known in the literature devoted to this 
subject (see e.g. Bensoussan et al. for a review on homogenization in periodic structures or 
Jikov et al. |13| for a monograph on stochastic homogenization). We will detail a few of them in 
the sequel of the paper. 

We finally mention that Theorem 13. II is very close to Theorem 3.1 in Guillotin-Plantard and Schott 
|lUj . Indeed, any dynamical (E, A, fi, T), as considered by the previous authors, generates an homo- 
geneous Markov chain with degenerate transitions: V/c S E, P(k, Tk) = 1 (but the state space E is 
then very large). In this framework, the condition J E fjdfi = {2d)^ 1 in Theorem 3.1 in jlUj implies 
that the expectation against \i of the transitions p(T l k, ej) = fj(T l x) and p{Tj:, —ej) = 1/d— fj(T l x) 
given in Subsection 3.1 in ^U] vanishes, and thus amounts in some sense to Assumption (A. 4) in our 
paper. In the same way, condition (H) in Theorem 3.1 in JOJ refers more or less to a "degenerate" 
central limit theorem for the underlying dynamical system and thus to our standing Assumption 
(A.3). 

However, our own Theorem 13 . 1 1 does not recover exactly the result in Guillotin-Plantard and Schott 
10 : due to the degeneracy of the transitions of a dynamical system, condition (A.3) cannot be 
satisfied. Actually, the reader must understand that condition (A.3) is a simple but very strong 
technical condition to ensure the validity of the so-called "central limit theorem" for Markov chains. 
It thus permits to draw up a clear framework in which the stochastic homogenization theory ap- 
plies, but is obviously far from being optimal (refer to Olla |22| for a complete overview on central 
limit theorems for Markov chains). 

Theorem 13. II is proved in Section 

3.2 Asymptotic Behaviour of the Reflected Walk 

The rescaled walk (Rf = m~^R^ 2t i )t>o satisfies the following "reflected version" of Theorem 

rm 

Theorem 3.2 The process (R[ m ^)t>o converges in law in B(]R + ,R d ) endowed with the Skorokhod 
topology towards the (unique) solution of the martingale problem with normal reflection on the 
boundary of the hypercube [0, l] d , with zero as initial condition and with C as underlying operator. 

The reflected SDE associated to the hypercube [0, \} d and to the operator C writes: 

dX t = b(X t )dt + dH t - dK t + a{X t )dB t . (3.2) 

We explain now the meaning of the different terms in the r.h.s of (|3.2j) . The quantity b(Xt)dt refers 
to the drift b(Xt)dt in (|3.1|) and a(Xt)dBt to the stochastic noise a(Xt)dBt in (|3.1jl . The new terms 
dH and dK stand for the differential elements of two continuous adapted processes with bounded 
variation that prevent X to leave the hypercube [0, l] d . More precisely, H l (the first coordinate of 
H) is an increasing process that pushes (when necessary) the process X 1 to the right to keep it 
above zero. In the same way, K 1 is an increasing process that pushes (again when necessary) the 
process X 1 to the left to keep it below one. For each £ S {2, . . . ,d}, H and K e act similarly in 
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the direction ee. Both processes H and K act in a minimal way. In particular, there is no reason 
to push the process X 1 when away from the boundary. 

The least action principle for H and K can be summarized as follows. The process iJ , for I G 
{1, . . . , d}, does not increase when X e is different from zero, and the process K , for £ € {1, . . . , d}, 
does not increase when X is different from one, i.e., 



For a complete review on reflected SDE's, we refer the reader to Tanaka [HJ, Lions and Sznitman 
[El and Saisho [2E]. 

The following corollary describes the asymptotic behaviour of the deadlock time and the deadlock 
point of the algorithm. We assume the matrix a — gg l to be uniformly elliptic to ensure a to be so. 
In short, the ellipticity of a avoids any singular behaviour (up to a P-null set) of the trajectories of 
the limit process X. 

Corollary 3.3 Assume that A G [1, d) and that the matrix a — gg l is uniformly elliptic. Denote 
by T the deadlock time of the limit process X: T = inf{t > 0, Xt £ Fq}, where Fq is given in 
Subsection \2.3L Then, the sequence of rescaled deadlock times m~ 2 T^ m ^ converges in law towards 
T. In the same way, the sequence of rescaled deadlock points m~ 1 R T ( m ) converges in law towards 



Theorem 13.21 and Corollary 13.31 are proved in Sectional 

3.3 Behaviour of the Limit Process in Dimension Two 

A crucial question for numerical applications consists in estimating precisely the mean value of 
the deadlock time of the limit system. When the matrix a is constant and diagonal and the drift 
b reduces to zero, several explicit computations in terms of Bessel functions are conceivable (see 
again [10]). 

In our more general framework, the story is quite different. In Delarue [3], we manage to establish 
in the two-dimensional case (i.e. d = 2) relevant estimates of the expectation of T (now denoted 
by T\ = inf{t > 0, X\ + Xf > A} to take into account the parameter A), and in particular 
to distinguish three different asymptotic regimes as the parameter A tends to two, each of these 
regimes depending on the covariance matrix a(l,l), and more precisely, on the sign of its off- 
diagonal components. 

Since the matrix a — gg l is assumed to be uniformly elliptic, the matrix a(l, 1) is positive and 
writes: 




(3.3) 



Xt . 




(3.4) 



with pi,p2 > and s e] — 1, 1[. The matrix a(l, 1) admits two eigenvalues: 




~ [pi + Pi + 5] , A 2 = \ [p\ +pj-S], where 5={p\ + p A 2 - 2(1 - 2s 2 )p 2 lP 2 2 ) 1/2 . (3.5) 
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Denote by E\ and E% the associated eigenvectors (up to a multiplicative constant). For s / 0, 

El= \{2sp xP2 )-\5 + pl-pf ) 

Since 5 + p\ — p\ > 0, the signs of the non-trivial coordinates of E\ and E% are given by the sign of 
s. The main eigenvector (i.e. E\) has two positive components for s > 0, and a positive one and a 
negative one for s < 0. Of course, if s vanishes, E\ and E<i reduce to the vectors of the canonical 
basis. 

The three different regimes can be distinguished as follows: 
Positive Case. 

If s > 0, the main eigenvector of a(l, 1) (i.e. E\) is globally oriented from to the neighbourhood 
of the corner (1,1) and tends to push the limit reflected diffusion towards the border line. The 
reflection on the boundary cancels most of the effects of the second eigenvalue and keeps on bringing 
back the diffusion along the main axis. As a consequence, the hitting time of the border line is 
rather small and the following asymptotic holds for the diffusion starting from 0: 

sup E(T A ) < +oo. 

1<A<2 

This phenomenon is illustrated below (see Figure when b reduces to and a is the constant 
matrix given by p\ = p 2 = 1 and s = 0,9. We have plotted there a simulated trajectory of the 
reflected diffusion process, starting from at time 0, and running from time to time 10 in the 
box [0, l] 2 . The algorithm used to simulate the reflected process is given in Slomihski [25]. The 
eigenvector E\ exactly matches (1,1)*. 



(2spip 2 ) 



+ P 2 2 



Pi 




Figure 3: One trajectory of the reflected process over [0, 10] for s = 0, 9. 



Negative Case. 

If s < 0, the main eigenvector of a(l,l) is globally oriented from (1,0) to the neighbourhood of 
the corner (0, 1) and attracts the diffusion away from the border line. Again, the reflection on the 
boundary cancels most of the effects of the second eigenvalue, and thus, acts now as a trap: the 
diffusion stays for a long time along the main axis and hardly touches the boundary. The hitting 
time satisfies the following asymptotic behaviour when the diffusion starts from 0: 

3ci, c 2 > 1, VA e]l, 2[, c^(2 - A)" C2 - a < E(T A ) < ci(2 - A)^ + c v 
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This point is illustrated by Figure below when b vanishes and a reduces to the constant matrix 
Pi = P2 = 1 an d s = —0,9 (again, the initial condition of the process is 0). The eigenvector E\ is 
given, in this case, by (1, —1)*. 




Figure 4: One trajectory over [0, 10] for s = —0, 9 (the trajectory is stopped before T). 



Null Case. 

The case s = is intermediate. Eigenvectors are parallel to the axes and the behaviour of the 
diffusion is close to the behaviour of the two-dimensional Brownian motion. For the initial condition 
0: 

3ci > 1, VA €]1, 2[, -q 1 ln(2 — A) — ei < E(T A ) < -cj ln(2 - A) + a. 

This is illustrated by Figure El below when b vanishes and a reduces to the identity matrix (the 
initial condition of the process is 0). 




Figure 5: One Trajectory over [0, 10] for s = 0. 



The following theorem sums up these different cases (we refer to Delarue 4 for the whole proof, 
and just provide a sketch of it in Sectional): 

Theorem 3.4 Assume that a — gg l is uniformly elliptic. Then, there exists a constant C j^^j > 1, 
depending only on known parameters and on the ellipticity constant of a, such that: 

1. For s > 0, sup E(7» < Qo 

le]l,2[ ] — tl 

2. For s < 0, set /?_ = -s > 0, 0+ = s(s - 3)(1 + s) _1 > 0. Then, 

W 2[, Qg<2 - £)-*- - < E(T t ) < Q^2 - + q^j 
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3. 7/ S = 0,WG]l,2[, -<^ln(2-^)-q^<E(r,)<-q^|ln(2-^) + C^| 

Note that Theorem 13.41 leaves open many questions. For example, we do not know how to compute, 
for s < 0, the exact value of the "true" exponent (5 = inf{c > 0, sup £g n 2 r [(2 — ^) C E(T))] < 
+00} . We are even unable to precise the asymptotic behaviour of f3 as s — * — 1 (note indeed that 
lim a _*_i/3_ = 1, lim s ^_i/? + = +00). 

We have very few ideas about the extension of Theorem 13.41 to the upper dimensional cases. The 
only accessible case for us is a(l, . . . , 1) = Id, Id denoting the identity matrix of size d: in this case, 
the analysis derives from the transience properties of the Brownian motion in dimension d > 3. 



4 Proof of Theorem 13.11 

The strategy of proof of Theorem 13.11 is well-known. We aim at writing the sequence m>l 
as a sequence of martingales, or at least of semimartingales with relevant drifts. The asymptotic 
behaviour of the martingale parts then follows from a classical central limit argument. 



4.1 Semimartingale Expansion of the Rescaled Walk 

(First Step. A First Martingale Form.) 

Denote by (J r n ) n >o the filtration (Jn' J )n>o and write (S n ) n >o as follows: 

S n +i = S n + J n +i 

= S n + E(j n+ i|J r n ) + J n+ i — E(J n+ i|jF n ). 

The increment Y® +1 = J n +i — E(J n+ i|^ r n ) appears as a martingale increment. Referring to (|2,2j) 
and to the definition of g given in Assumption (A. 4), the underlying conditional expectation writes 

E(j n+1 |^ n ) =g(£ n ,S n /m). (4.1) 

Thus, S n+ i now writes: 

S n+1 = S n + g(£ n , S n /m) + y n ° +1 . (4.2) 
We deduce the following expression for the rescaled walk: 

[m 2 t]-l [m 2 t] 

St^m- 1 £ g[^S k /m)+m- l Y, Y k- (4-3) 

k=0 k=l 

Note at this early stage of the proof that condition (A. 4) is necessary to establish Theorem 13.11 
Assume indeed that g does not depend on S n /m. Then, the first term in the above right hand side 
reduces to 

\m 2 t\-l 

m ^ E ( 4 - 4 ) 

fc=0 

The ergodic theorem for Markov chains then yields: 

\m 2 t\-l 

mT 2 \^ Q{£k) —> t I qdri as m —> +00. 



fc=o Je 
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In particular, if the expectation of g with respect to \i doesn't vanish, the term Q4.4|) explodes with 
m. The question is then to study the asymptotic behaviour of (|4,4j) under the centering condition 
(A. 4), or in other words to establish a central limit theorem for the sequence (g(^ n ))n>o- 

(Second Step. Auxiliary Problems.) 

The basic strategy consists in solving the system of Poisson equations driven by the generator of 
the chain (£n)n>o an d by the M. d valued drift g. In our framework, this system extends to a family 
of systems since g does depend on the extra parameter y G M. d . We thus investigate: 

(I-P)v = g(-,y), y£R d . (4.5) 

For y G M. d , the solvability of this equation is ensured by the Doeblin condition (A. 3). Referring to 
Remark 2, Theorem 4.3.18 in Dacunha-Castelle and Duflo [3], we derive from (A. 3) that there exist 
two constants T and c, < c < 1, depending only on known parameters appearing in Assumptions 
(A. 1-4), such that: 

Vn > 0, Vi G E, \P n (iJ)-Kj)\ < ^ n . (4.6) 
We then claim (see e.g. Theorem 4.3.18 in Dacunha-Castelle and Duflo 3 ): 

Proposition 4.1 Under (A. 1-4), the equation j4-5\ ) is solvable for every y G Moreover, there 
is a unique centered solution with respect to the measure fx. It is given by: 

Wi G E, v(i, y) = J2 E [ pn &Mj, V)] ■ 

n>0j£E 

Note carefully that v (i, y) belongs to M. d . Note also from (A. 4) and (J4.6|) that v(i, y) is well defined 
since 

v{i,y) = Y J Y,[ pn ^^-^Mj,y)- 

n>0j£E 

The question of the regularity of g and v turns out to be crucial in the statement of Theorem 13.11 
The following Lemmas directly follow from Assumption (A.l) and (|4.6j) : 

Lemma 4.2 The function g(i, y) is bounded by 1 and is twice continuously differentiable with 
respect to y with bounded derivatives. In particular, g is Lipschitz continuous in y, uniformly in i. 
The Lipschitz constant is denoted by Cj^jJ 

Lemma 4.3 The function v(i,y) is bounded by T(l — c) _1 and is twice continuously differentiable 
in y with bounded derivatives. In particular, v is Lipschitz continuous in y, uniformly in i. The 
Lipschitz constant is denoted by 

(Third Step. Final Semimartingale Form.) 
We claim: 

Proposition 4.4 Define for (i, y) G E x IR^: 

b {m) (i,y) = ^2[p(i,y,u)((v - g)(i,y + u/m) - (v - g)(i, y))] . 
uev 
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and denote by (Y n ) n >i the sequence of corrected martingale increments: 

Y n = J n + v(£ n ,S n /m) - E[j n + v(£ n , Sn/rnjlfn-i]. 

Then, the walk (S n ) n >o writes: 

S n+ i + v(S n+1 ,S n +i/m) =S n + v(£ n , S n /m) + 6 (m) (£ n , S n /m) + Y n+1 . (4.7) 

In particular, the process (S n + v(£ n , S n /m)) n >o is a semimartingale. The term b^ appears as a 
drift increment. It satisfies: Vm > l,i 6 E, y G M. d , \b^ m \i,y)\ < C ^j^ m, where C ^j^j depends 
only on known parameters in (A. 1-4). 

Proof. Apply v to the couple (£ n +lj Sn+i/m) f° r a given n > 0. In this perspective, note that: 
E(v(£ n+ x,S n+ i/m)\F n ) = ^2 [ p (Cn,i)p(U, S n /m,u)v(i, (S n + u)/m)] 

i£E,ueV 

= ^2 p{in,S n /m,u)[ ^P(i n ,i)v{i, {S n + u) I W) 



uev 



According to Proposition 14.11 we have: 

^2[P{Cn,i)v(i,{S n + u)/m)] = (v-g)(£ n , (S n + u)/m), 

yielding: 
E(v (^+i, Sn+x/m)]^) 

= (V ~ g)(Cn, S n /m) + ^ 



p{in,S n /m,u){ (v - g)[£ n , (S n + u)/m) - (v - g)(£, n ,S n /m} 



Due to the definition of b^ m \ this writes in an equivalent form: 

g(U,S n /m) =v(Cn,S n /m) -E(v(£ n+1 ,S n+1 /m)\r n ) +b^^ m S n /m). (4.8) 
From (|4~^|) and (|4~H|) . derive that: 

S n +1 = S n + Jn+1 + v(^n, S n /m) — E(j ra+ i + u(£ n +l) 

Sn + l/m)\^n)+b {m) {in,Sn/m). 

Thanks to the definition of (Y n ) n >i, we recover the required writing. The bound for &( m ) derives 
from Lemmas 14.21 and 14.31 (Xipschitz continuity of g and v). □ 

Investigate the bracket of the martingale part of (S n + v(£ n , S n /m)) n >o- 

Proposition 4.5 Denote by (M n ) n >o the square integrable martingale: M n = 5^fc=i^fc- Then, the 
bracket of (M n ) n >i is given by: 

n 

(M)n = Y, aim) (Zk-i,S k ^/m), 
k=l 
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with (see the statement of Theorem \°J.l\ for the definition of a): 
(£„, S n /m) = [a-(v + b^)(v + b^f] (£ n , S n /m) 



(4.9) 



Moreover, there exists a finite constant depending only on known parameters in (A. 1-4), 

such that: Vm > 1, i E E, y E R d , \a^(i,y)\ < 

The notation (Pv){i,y) in Q4.9|) (the analogue holds for (P(vv t ))(i,y)) stands for the <i-dimensional 
vector (Pvi(i,y), . . . ,Pv d (i,y)y. 

Proof. Recall that \J n \ is almost surely bounded by 1 and that v is a bounded function (see Lemma 
14.31) . Hence, the variables (Y n ) n >i are almost surely bounded and (M n ) n >i is square-integrable. 
The bracket is given by: 

n 

(M) n = ^[Jk + v{&, Sfc/m)|^ fc _i] , 
fc=l 

where V[.|.Ffc_i] denotes the conditional covariance with respect to the cr-field lFk—i- Due to (|4.1() 
and (|4.8|) . we have: 

®[Jn+l+v{£n+l,S n+ l/m)\T n ] = (v + b^)(C n ,S n /m). 

Since V[Z\F n ] = E[ZZ*|^ n ] - E[Z|.F n ]E[Z*|.F n ], 
V[J n+ i +t;(£ n+ i,S' n+ i/m)|j n ] 

= E ( J n+1 + v(£ n +i, S n+1 /m)) (J n+1 + v(f n+1 , S , n+ i/m)) i |^ n 

- [( V + 6 {m) )(^ + 6 (m) )*](^,5 n /m) 

= P(^,i)jp(^ n ,5 n /m,u)[it + u(i, (S n + u)/m)] [u + v(i, (S n + u)/m)] 1 

i€E ngV L 

- [(v + b^)(v + b^) t ](Cn,S n /m). 
Reducing the sum over i E E: 

P^n,S n /m,U)UU t ] 

uev 



uev 



uev 

+ H^ n ' 5 ' n / m ' u ) u ( Pv ) t {^.n, (S n + u)/m)] 

P(€n, S n /m,u)(Pv)(£ n , (S n + u)/m)u t 

P(£n, Sn/m, u) [P(W)] (Cn, (£ n + u)/m 



- [(^ + 6 (m) )(t; + 6 (m) )*](^,5 n /m). 
This recovers the form of a^ m \ The bound for follows from Lemmas I4.2I and I4.3I □ 
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4.2 Scaling Procedure 

(First Step. Semimartingale Form of the Rescaled Process.) 
We deduce from Proposition 14.41 that for all t > : 

\m 2 t\-l im 2 t\ 

5| m) +m-% 2iJ ,5< ra) )=m- 1 ^,0)+ m ^ 1 £ 6<"0(&, S k /m) + m' 1 ^ Y fc . 

fc=0 fc=l 

Of course, if |jn 2 £j = 0, the latter sums reduce to zero. 

Since the function v is bounded (see Lemma l4.3|) . the supremum norm of the process {'m~ 1 v(^ m 2 t ^ , 

St m ^))t>o tends almost surely towards as m increases. In particular, writes, up to a negligible 
term: 

\m 2 t\-l \m 2 t\ 

Vt>0, Sl m) =m- 1 ]T ^ m \i k ,S k /m)+m- 1 £ Y k + 0(l/m). (4.10) 

k=0 k=l 

Set now for the sake of simplicity: 

\m 2 t\-l \m?t\ 

B[ m) ^ m- 1 £ b^(C k ,S k /m), = m' 1 £ Y k . (4.11) 

A:=0 fe=l 

(Second Step. Tightness Properties.) 

Proposition 4.6 The family (S (m \ M^, (M^)) m >i is C-tight in B(R+, R 2d+d2 ). 

Proof. We first establish the tightness of the processes ) m >i i n the space of cad-lag 

functions endowed with the Skorokhod topology. To this end, note from Proposition 14.41 (estimates 
for (&( m )) m >i) that for any (s,t) G (M+) 2 , < s < t: 

lB (m) _ < Cm -2 {[m 2 t] _ [m 2 s]) 

<C((t-s) + m- 2 ). 

Therefore, the family (-B^ m - ) ) m >i is C-tight in D(M + ,lR a! ) endowed with the Skorokhod topology: all 
the limit processes of the sequence have continuous trajectories (see e.g. Jacod and Shiryaev [T2*] . 
Propositions VI 3.25 and VI 3.26). 

We turn now to the process . It is a square integrable martingale with respect to the filtration 
(J~\ m 2 t\ )t>o- Moreover, from Proposition 14. 5[ the bracket of the martingale is given by: 

[m 2 t\-l 

(M^) t = m- 2 £ a^(^ k ,S k /m). (4.13) 

k=0 

Since the functions (a( m )) m >i are uniformly bounded, the family ((M( m ))) m >i is also C-tight in 
, R M ), by the same argument as above. 



Referring to Theorem VI 4.13 in Jacod and Shiryaev ^2], we deduce that the family (M( m -*) m >i 
is tight in B(M + ,]R d ). It is even C-tight since the jumps of the martingale M( m ) are bounded in 
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absolute value by C/m. 

Finally, we conclude that the family (S^, M< m ), (M( m ))) m >i is C-tight in 
Corollary VI 3.33 in Jacod and Shiryaev |12]). □ 

(Third Step. Extraction of a Converging Subsequence.) 



(see 



Lemma 4.7 Every limit (X,M,L) (in the sense of the weak convergence in D(R + , M 2a!+d2 ) y ) of 
the sequence of processes (S*( m ),M( m ),(M( m ))) m >i is continuous. Moreover, M is a square inte- 
grable martingale with respect to the right continuous filtration generated by the process (X,M,L) 
augmented with F-null sets. Its bracket coincides with L. 



Proof. According to Proposition IX 1.17 in Jacod and Shiryaev [12], M is a square integrable 
martingale. The same argument shows that M 2 — L is a martingale. It is readily seen from (|4.13f) 
that L is of bounded variation and satisfies Lq = 0. Thus, (M) = L (see Theorem I 4.2 in |12j). □ 

We aim at expressing the process (M). To this end, we have to study the asymptotic behaviour of 

(IW): 

Lemma 4.8 Referring to the definition of a in the statement of Theorem \H.l\ 

\m 2 t\-l 



sup 

0<t<T 



m)\ 



m 



E a(S k /m) 



k=0 



-probability 



0asm-> +oo. 



Proof. The function from (|4.9j) uniformly converges on E x M d towards the function a given 
by: 



Vi 6 E, Vy € R d , a(i, y) = a(i, y) + [giJPvf + (Pv)g t + P(«u*) - vv l ] (i, y) 

= a(i, y) + [gv* + vg* - 2gg t + P(uv f ) - w*] (i, y). 



(4.14) 



This follows from the a priori estimates for (M m **) m >i (see Proposition 14. 4|) and from Lemma 14.31 
Thus, (M( m )) behaves like: 

\m 2 t\-l 

(m- 2 Y a fa,S k /m)) . (4.15) 
fc=o ' * >0 



Lemma 14.81 then follows from Proposition 14.91 below. □ 

Proposition 4.9 Let f : E x W 1 — > R 6e bounded and uniformly Lipschitz continuous with respect 
to the second variable. Then, we have for all T > 0: 



sup 

o<t<r 



m 



[m 2 tj-l 



f(£k,Sk/m) 



k=0 



\m z t\~l 
fc=0 



-probability 



0asm-> +oo, 



with Vy G R d , /(y) 



f(i,y)dfj,(i). 



The proof of Proposition 14.91 relies on the ergodic theorem for (£ n )n>o an d ° n the tightness of the 
family (S^) 

m >i- It is postponed to the end of the section. 
(Third Step. Identification of the Limit.) 
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Lemma 4.10 Under the assumptions and notations of Lemmas \4- 7| and \4-8\ and of Proposition 
\4-9\ the bracket of M is given by Vt > 0, {M)t = Jq a(X s )ds, and the couple (X,M) satisfies: 

Vt>0, X t = f b(X s )ds + M t - (4.16) 
Jo 

Proof. Focus first on the limit form of the bracket. Note to this end that for every t > 0: 

[m 2 t]-l 



m 2 



V d(S k /m)= [ a(S^)ds + 0(l/m). (4.17) 
k=o Jo 



The convergence will then follow from the one of (S^ m ') m >o. 

In fact, Proposition 14.91 applies to studying the asymptotic behaviour of the drift part. Derive 
indeed from Lemmas 14.21 and 14.31 and Proposition 14.41 that mb( m ' uniformly converges on E x M. d 
towards (V y v — V y g)g. Thus, the asymptotic behaviour of in (|4.11j) just reduces to the one 
of: 

L^J-i d d , v i_ gi ^ 
m 2 Yl Htk,S k /m), b = (V y v-V y g)g=(Y ^ 9j) l<i<d - 



k=0 j=l 

According once again to Proposition 14.91 this just amounts to study: 



\m 2 t\-l 

-2 



m ' Y KSk/m)= [ b(Si m) )ds + 0(l/m), 

k=o Jo (4.18) 



b (y)= / b (hy) d K i ) = / [(V y f - V y g)g](i,y)d(j,(i). 

JE JE 

Note now from Lemma 14.71 and from the continuity of a and b (see Lemma 14.31 and the proof of 
Lemma H~H)l that the processes: 

5 (m) ,M (m) , [ b(Si m) )ds, [ d(Si m) )ds] (4.19) 

JO JO J m>l 

converge in law in D(M+, R 3d+d2 ) towards (x,M, J b(X s )ds, J d(X s )ds^j. 

From Lemmas H~7I and Ol and (gTTTJ), the bracket of M is given by: Vt > 0, (M) t = f t d(X s )ds, 
and the relation (|4~TU|) yields (jUTHJ). □ 

We are now in position to complete the proof of Theorem l3.ll Thanks to Lemmas 14. 21 and 14.31 the 
coefficients a and b given in the statement of Theorem 13.11 are once differentiable with bounded 
derivatives. It is then well known that X in Lemma 14. 1UI turns to be the unique solution (starting 
from the origin at time zero) to the martingale problem associated to (b, a) (see e.g. Chapter VI 
in Stroock and Varadhan |3Uj ) . 

We now prove that the diffusion matrix a is elliptic when a — gg l is uniformly elliptic. Referring 
to the definition of a (see Theorem 13.1(1 . it is sufficient to prove that the symmetric matrix: 

[vg* + gv* - gg t ](i,y)dfi(i), 
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is non-negative. Recall that g = v — Pv. Thus: 

vg* + gv l - gg l = vv* - v(Pv) t + vv l - {Pv)v l - vv l - (Pv)(Pv) t + v(Pv) t + (Pv)u* 
= vv t - (Pv)(Pv) t . 

We see from Schwarz inequality that for every (x,y) G R d x R d : 

(x,(Pv)(Pv) t x)(i,y)df,(i)= [ (x,Pv) 2 (i,y)d^i) 

E JE 



(4.20) 



< I ^2P(i,j)(x,v(j,y)) 2 d^i) (4.21) 

•J E Tp 



j£E 

(x,v(i,y)} 2 d^(i) 

E 

(x, (vv t )x)(i,y)dn(i). 

E 

Thanks to (|4.2Uj) and (|4.21|) , we complete the proof of Theorem 13.11 □ 
4.3 Proof of Proposition 14.91 

The proof follows the strategy in Pardoux and Veretennikov jH]. For this reason, we just present 
a sketch of it. 

Let T be an arbitrary positive real number. Recall that (<S( m )) m >i is tight in D(R + ,M d ). Hence, 
for a given e > 0, there exists a compact set A C ID(IR+,IR rf ) such that, for every m > 1, 
belongs to A up to the probability e: 

P{S (m) G A} > 1 - s. (4.22) 

Since A is compact for the Skorokhod topology, we know, for 5 small enough, that every function 
£ £ A satisfies: 

w' T+1 (£,6) <e, (4.23) 

where w' T , x denotes the usual modulus of continuity on [0, T + 1] of a cad-lag function (see e.g. 
Jacod and Shiryaev, VI 1.8 and Theorem VI 1.14). Moreover, there exist £%, ■ ■ ■ ,Inq i n A such 
that, for every i G A, we can find an integer i G {1, . . . ,iVo} and a strictly increasing continuous 
function Aj from [0, T] into [0, T + 1/2] such that: 

sup [\t - Xi(t)\] < 5 2 , sup [\£(t) - £i(Xi(t)) |] < 5 2 . (4.24) 
te[o,T] te[o,T] 

Since ii satisfies ()4.23(l . there exists a step function yi : [0, T + 1/2] — > R, with all steps but maybe 
the last one of length greater than 5, such that the supremum norm between y» and ii is less than 
e on [0,T + 1/2]. The second term in P^lj) then yields: sup tg[0jT ] [\£(t) - yi(\i(t))\] < 5 2 + e. 
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Choose now u G Q such that S^ m '(u}) G A. Deduce that there exists i G {1, . . . , ./Vo} such that the 
supremum on [0, T] of the distance between £( m ) and ?/j o Aj is less than e + 5 2 . In particular: 

sup Im-^fc-^AiCfc/m 2 ))! < sup |5 t (m) - W (Ai(t))| < e + 8 2 . 

fce[0,|m 2 T-lJ] t£[0,T] 

Since / is Lipschitz continuous, 

Lm 2 tJ-l 



sup 

0<t<T 



m 



E |/&><Vm) - /(£fc,y*(Wm 2 )))| 



fc=0 



< CT(e + <5 2 ). 



(4.25) 



Define now by to = 0, t\, . . . , t p = T a subdivision of [0, T] with respect to the steps of j/j. Recall 
that - > 5 for j G {1, ... ,p - 1}. According to (jOijl . for f G G = U^[t fc _i + <5 2 ,i fc - 
(5 2 ]U]t p _i+5 2 AT, T], Aj(i) and t belong to the same class of the subdivision, so that yi(Aj(i)) = Vi{t). 
Hence, 



lim sup 



[m 2 T\-l 



E \f(Zk,yi{\(k/m 2 ))) -/(&,w(fc/ro 2 ))| 



k=0 



< C lim sup 

m— >+oo 



Lm 2 TJ-l 



-2 



E ^t^gOA" 2 



fc=0 



(4.26) 



C(T - //(G)) < 2CS 2 p < 2CT5. 



It follows from (|4.22|) . (|4.25|) and Ij4.26|) that there exists a constant C (not depending on <5, e, m 
and iVg) such that for m large enough: 



U 

ie{i,...,JV } 



sup 

0<t<T 



\m 2 t\-l 



m 



E |/(&>&/™)-/(&^i(V™ 2 ))| 



fc=0 



> C(<y + e) M < e. (4.27) 



Note that the same estimate holds with Sk/m) — f{£,k-,yi{k/m 2 ))\ replaced by |/(£fc, Sk/m) 

f(^,yi(k/m 2 ))\ + \f{S k /m) - f(y t (k/m 2 ))\. Moreover, for all % G Ex 

[m 2 T] p |m 2 %J 

E ffaMk/m 2 )) =J2 E f{tk, yi {k/m 2 )) 

3=1 fc=|m%_ij + l 



k=l 



[m 2 tj 



E E HtkMtj-i))- 

3=1 fc=|m 2 t,_ij + l 



The ergodic theorem for (£ n )n>o yields: 

\m 2 T\ 



m 



E fiikMk/™ 2 )) ^ E[fe - = / Ryi{s))ds. 

k=i 3=1 Jo 
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Obviously, m 2 ^i=o^ 1 f{yi{k/m 2 ) ~^ Jq f(Ui( s ))ds as m tends to +00. Hence, for m large 
enough: 



\m 2 T\-l 



[m 2 TJ-l 



m 



m 



Y f(m(k/m 2 )) 



k=0 



k=0 



This actually holds for every t £ [0,T]. Choosing t of the form t = ke, k £ {0, . . . , [e 1 TJ}, it 
comes (for m large): 



sup 

te{o,...,£Le- 1 rj} 



[m 2 t\-l 



[m 2 t\-l 



m 



E f^kMk/m 2 )) 



m 



Kviik/m 2 )) 



k=0 k=0 

Since / is bounded, it comes finally for a suitable constant C > (and for m large): 



D < sup 
[o<t<T 



\_m 2 t\-l 



[m 2 t\-l 



in 



-2 



Y fttkMk/m 2 )) 



m 



Y Km(k/m 2 )) 



k=0 



k=0 



> Ce \ < TNq s, (4.28) 



The sum over i £ {1, . . . , No} in (|4.28|) is bounded by e. Combining (|4.27|) and (|4.28j) . we complete 
the proof. □ 



5 Proofs of Theorem 13.21 and Corollary 13.31 
5.1 Proof of Theorem GO 

The proof of Theorem 13.21 is rather similar to the one of Theorem 13.11 The main difference comes 
from the definition of the mean function g. In the new setting, it writes for (i, y) £ E x [0, l] d 
h(i, y) = YlueV u< l(^ Vi u )- ^ n particular, for £ £ {1, . . . , d}: 



hih y) = 9e(i, y), for <y £ <l, 

he(i,y) =p(i,y,e e ) + p(i,y,-ee), for y e = 0, 

ht{i,y) = -(p(i,y,ee) + p(i,y,-e e )), for y e = 1, 



(5.1) 



where gi (resp. hi) denotes the £th coordinate of g (resp. h). The process R satisfies the analogue 
of (|4,3|) with g replaced by h. Due to (|5.1|) . the drift hi, for £ £ {1, . . . , d}, pushes upwards the £th 
coordinate of R when matching 0, and downwards when matching 1. 

The drifts of the processes R and S just differ in the extra term (h — g). The following lemma 
derives from (|5.1|) : 

Lemma 5.1 For £ £ {1, . . . , d} and y £ [0, l] d , ye = =>- (h — g)i > 0, yt = 1 =>• (h — g)i < and 
y t e]0,l[=*(h-g) t = 0. 

(First Step. Semimartingale Form.) 

Of course, the function h does not satisfy the centering condition (A. 4). However, for g as in 
Section 01 we can still consider the solution v to the family of Poisson equations (|4.5|) . Following 
the proof of Proposition 14.41 it is plain to derive: 
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Proposition 5.2 Define as in Proposition \4-4[ but with p(i,y,u) replaced by q(i,y,u). 
Then, the walk (R n ) n >o satisfies: 

R n +i + v{£ n+ i,R n+ i/m) = R n + v(£ n ,R n /m) + (h-g)(£ n ,R n /m) 

+ c ( ~ m \U,Rn/m) + Z n+u 

where E(Z n +i| J-^' R ) = 0. The drift increment c^ satisfies the same bound as b^ m h Vm > l,i 6 
E,yeR d , \c^(i,y)\ <0£gm. 

We are then able to control the variation of R: 

Lemma 5.3 Define e = e± + - ■ - + ed- Then, there exists a constant depending only on known 

parameters in (A. 1-4), such that for 1 < n < p < q, 



(5.3) 



(5.4) 



q-l q-l 

\Rg - R n \ 2 + 2^(me - R n , (h - g)~^ k , R k /m)) + 2 ^(i? n , (h - g) + (^ k , R k /m)) 

k=p k=p 

9-1 

< \Rp - R n \ 2 + Q^jg -p + m) + 2^(R k - R n ,Z k+1 ). 

k=p 

Proof. For k £ {n, . . . , q — 1}, derive from ()5.2|) : 

\Rk+l — Rn\ 2 = \Rk — Rn\ 2 + 2(-Rfc+l — Rk, Rk ~ Rn) + \Rk+l ~ Rk\ 2 

= \R k - R n \ 2 + 2(R k - R^, (h - g)(£ k ,Rk/m)) 
- 2(R k - R n ,v{^ k+1 ,R k+1 /m) - v(£ k ,R k /m)) 
+ 2(R k - R n , c^ Rk/m)) + 2(R k -R n ,Z k+1 ) 
+ \Rk+l — Rk\ 2 - 

= \R k - R n \ 2 + 2T(1, k) + 2T(2, k) + 2T(3, k) + 2T(4, k) + T(5, k). 

Apply Lemma 15.11 and derive that —T(l,k) writes (me — R n ,(h — g)~ (£ k , R k /m)) + (R n ,(h — 
9) + (^k,Rk/nT'))- Summing T(l, k), for k running from p to q — 1, we obtain the second and third 
terms in the l.h.s of (|5.3j) . 

Perform now an Abel transform to handle ^ q }Z}pT(2,k). The boundary conditions are bounded 
by Cm since R takes its values in [0, m] d and i; is bounded (see Lemma I4.3JI . Moreover, the sum 
Yl k =!p+i(Rk — Rk-i, v(^ k ,R k /m)) is bounded by C(q — p) since the jumps of R are bounded by 1. 
Hence, the term Yl k =pT(2, k) is bounded by C(q — p + m). 

Due to the bound for in Proposition 15. 2| the sum 2~^fc=p^(3 ; k) is bounded by C(q — p). The 
sum of T(4, /c) over A: provides the third term in the r.h.s of (|5.3|) . 

Finally, since the jumps of i? are bounded by 1, the sum of T(5, k) over A; is bounded by (q — p). 
□ 

We turn now to the martingale part in Q5.2[) . Following the proof of Proposition 14.51 we claim: 
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Proposition 5.4 Define the martingale N by Vn > 1, N n = 5^fc=l ^ s bracket {N) n satisfies 
a suitable version of Proposition \4-5\ with replaced by , given by: 



i (m) (U, Rn/m) =[a-(v + h-g + c^)(v + h-g + c^f] (£ n , R n /m) 

+ S "/ m ' u ) [ p ( vvt ) + u ( Pv f + (Pv)u l ] (£„, (#n + «)/m) 



(Second Step. Tightness Properties.) 

The scaling procedure then applies as follows. Define as in (|4.11jl : 



|m 2 ij-l 



£r( m ) 

f/ t = m 



k=0 
\m 2 t\-l 



1 ^ (h-g)+(t k ,R*/m), K 



(m) 



|_m 2 i]-l 



■/?) 



- 1 ^ (h-g)-(Z k ,R k /m) 



k=0 



k=0 



(5.5) 



Proposition 5.5 The family (R^ m \N^ m \ (N^), H^ m \ if( m )) m >i is C-%M in D(R+, M 4d+d2 ). 

Proof. Following Proposition EH the family (iV*( m ), (iV( m ))) m >i is C-tight in B(K+, R d+(f ). 

Apply now Lemma 15.31 to establish the same property for (^ m ')»n>i- Choose n = p = [m 2 s\ and 
q = [m 2 t\ , with s < t. Since all the terms in the l.h.s of ()5.3j) are nonnegative, it comes: 

\m 2 t\-l 

\RT ] -R^\ 2 <C(t-s) + 2m- 2 ^2 (Rk,Z k+1 ) + 0(l/m). (5.6) 

k= [m 2 sj 

The second term in the r.h.s of the above inequality defines a family of martingales. It is plain 
to see that it is C-tight in 1D)(M+,M) and thus to derive the C-tightness of the family (R^ m ') m >i in 
B(R+,R d ). 

We turn finally to the tightness of (H^ m \ K^) m >i. Choose p and q as above and n = in ()5.3|) . 
i.e. R n = 0. It comes: 

(e,Ri m) - RW) < -(i? f (m) ) 2 + [R^ f 

[m 2 t\-l 

+ C{t-s + l/m) + 2m- 2 ^ (R k ,Z k+1 ). 

k= \m 2 s\ 

Note that |(-Rj m) ) 2 - (i?i m) ) 2 | < 2|i?J m) - i?i m) | since R^ belongs to [0, l] d . Due to deduce 
that (i?( m )) m >i is C-tight in D(M + ,lR d ). Of course, the same holds for (#( m )) m >i. 

(Third Step. Extraction of a Converging Subsequence and Identification of the Limit.) 

Proposition 5.6 Any weak limit (X,N,H,K) of the sequence (R( m \ N^ m \ H^ m \ K^) m >i sat- 
isfies: 

Vt > 0, X t = I b(X s )ds + H t - K t + N t , (5.7) 
Jo 
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where N is a square integrable continuous martingale whose bracket writes (N)t = Jq a,(X s )ds, H 
and K are two nondecreasing continuous processes, matching at zero and satisfying condition 

Proof. Note first that Proposition 14.91 still applies in the reflected setting, with S replaced by R. 

Focus now on the asymptotic behaviour of the drift c^ in (|5.2[) . As mb^ does, mc( m ) uniformly 
converges towards (V y v — 'V y g)g. Following Q4.18JI . this provides the form of the limit drift. 

Concerning the martingale part, we now prove the analogue of Lemma 14.81 The point is to study 
the asymptotic form of the bracket (N^). It well seen from (|5.5j) that uniformly converges 
towards a (as in (|4.14|0 plus a corrector term that writes (g — h) r Y t + r y(g—h) for a suitable bounded 
function 7. Of course, the part in a satisfies Proposition 14.91 The second part vanishes since: 



\m 2 t\-l 

- £ \h-gm,R k /m)<m-'(Ht' , +K^" 
k=0 



This provides the limit form of the brackets ({N^ m ^)) m >i and thus the form of the limit bracket. 

We turn finally to the processes H and K. It is clear that Hq = Kq = and that H and K are 
nondecreasing. It thus remains to verify ()3.3j) . Focus on the case of H. Since (R^ m \ H^ m ^) m >i 
converges weakly (up to a subsequence) in the Polish space 10(R, R 2d ) towards (X,H), we can 
assume without loss of generality from the Skorokhod representation theorem (see Theorem 11.86. 1 
in Rogers and Williams 26 ) that the convergence holds almost surely. In particular, for T > 0, 
for £ G {1, ... ,d}, for almost every uj G $7, (H^ m \e£), seen as a distribution over [0, T] converges 
weakly towards {H,ee) seen as a distribution over [0, T]. Hence, for a continuous function / on 

[ °' T]: 

lim F f(t)d(H^\e e ) t = f f{t)dH l t . (5.8) 

Since X is continuous, Q5.8JI holds with f(Xt) instead of f(t). Recall moreover from Proposition 
VI. 1.17 in Jacod and Shiryaev J2] that sup < t <^ ~ Xt\ tends to almost surely. Since the 



supremum sup m>1 is almost surely finite, Q5.8|) yields: 

lim f T f(R?)d(S( m \e e ) t = ^ f{X t )dHl (5.9) 
n ^+°°Jo Jo 

From Lemma f5.ll for a continuous function / with support included in ]0, 1], the term f(Xf)dHf 
vanishes. It is plain to derive the first equality in 1)3. 3 J) . □ 

According to Proposition 15. 4[ (X, H, K) satisfies the martingale problem associated to (b, a) with 
normal reflection on d[0, l] d . Thanks to Theorem 3.1 in Lions and Sznitman [16] ■ this martingale 
problem is uniquely solvable. This completes the proof of Theorem 13.21 □ 

5.2 Proof of Corollary EP 

The proof of Corollary 13 .31 relies on the mapping Theorem 2.7, Chapter I, in Billingsley Indeed, 
consider the mapping ^ : x G B(]R+, M. d ) 1— > inf{t > 0,xj G Fq} (inf = +00). If x denotes a 
continuous function from M + into M. d and (x n ) n >\ a sequence of cad-lag functions from M + into W 1 
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converging towards x for the Skorokhod topology, we know from Proposition VI. 1.17 in Jacod and 
Shiryaev |12| that (ar n ) n >i converges uniformly towards x on compact subsets of R+. Assume now 
that ty(x) is finite and that for every rj > there exists t £]^(x), ^(x) + 77 [ such that x(t) belongs 
to the interior of Fq (the point x(^(x)) is then said to be regular). Then, ^(x n ) tends to 9(x). 
In particular, if almost every trajectory of X satisfies these conditions, ^{R( m >) converges in law 
towards W(X). 

To prove that almost every trajectory of X satisfies the latter conditions, consider the solution X to 
the martingale problem associated to the couple (b, a), where b and a are periodic functions of period 
two in each direction of the space given by (b(y),a(y)) = (N (y)b(H(y)) , a(II(y))) for y G [0,2[ rf , 
with U(y) = (vr(yi), . . . ,7r(y d )), tt(z) = zl [0<1] (z) + (2 - z)l] 1}2 ](z) and N(y) = diag(sign(l - 
y±), . . . ,sign(l — yd)) (diag(x) denotes the diagonal matrix whose diagonal is given by x). Due to 
the boundedness of b and to the continuity and to the ellipticity of a, such a martingale problem is 
uniquely solvable (see Stroock and Varadhan |30j . Chapter VII). Extend II from a 2Z d -periodicity 
argument to the whole set R d , and derive from the Ito-Tanaka formula that II(X) satisfies the 
martingale problem associated to (6, a) with normal reflection on d[0, l] d . In particular, H(X) and 
X have the same law. It is then sufficient to prove that almost every trajectory of 11(A) satisfies 
the conditions given in the above paragraph, and thus to prove that almost every trajectory of X 
hits F = U fee2 z< i (^ + -^b) an d that every point of the boundary of F is regular for X. 

Due to the boundedness of b and a and to the uniform ellipticity of a, deduce from Pardoux 
23 , Section 2, that the process X, seen as a diffusion process with values in the torus R a! /(2Z o! ) 
is recurrent, and from Pinsky Theorem 3.3, Section 3, Chapter II, that every point of the 
boundary of F is regular for X. 

The second assertion in Corollary 13.31 follows again from the mapping theorem. □ 

6 Sketch of the Proof of Theorem 13.41 

We now present several ideas to establish Theorem 13.41 Again, the whole proof is given in Delarue 
®- 

We denote by (A) the assumption satisfied by a and b under the statement of Theorem 13. 41 Briefly, 
a and b are bounded and Lipschitz continuous, and a is uniformly elliptic. 

6.1 Description of the Method 

As a starting point of the proof, recall that the two-dimensional Brownian motion B never hits zero 
at a positive time, but hits infinitely often any neighbourhood of zero with probability one. The 
proof of this result (see e.g. Friedman |5j) relies on the differential form of the Bessel process of index 
1, i.e. of the process \B\. In short, for B different from zero, d\B\ writes d\Bt\ = \/{2\Bt \)dt + dBt, 
where B denotes a one-dimensional Brownian motion. 

The common strategy to investigate the recurrence and transience properties of \B\ then consists in 
exhibiting a Lyapunov function for the process | B \ (see again Friedman [Sj for a complete review on 
this topic). In dimension two, i.e. in our specific setting, the function In is harmonic for the process 
\B\ (the Ito formula yields for B different from zero: dln(|St|) = \Bt\~ l dBt). In particular, for 
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|£>o| = 1) f° r a large N > 1 and for the stopping time tn = t n -iAt2, with t x = inf{i > 0, \Bt\ = x}: 

= E(ln(|£ o |))=E(m(|£ fjv |)) 

= - ln(JV)P{r|v-i < r 2 } + ln(2)P{r 2 < t n -i}. 

Hence, limjy^ +00 F{t n -i < T2} = 0, so that P{tq < T2} = 0. Of course, for every r > 2, the 
same property holds: IP{ro < r r } = 0. Letting r tend towards +00, it comes P{t < = 0, as 

announced above. 

More generally, if Y denotes a Bessel process of index 7 > -1, 7 / 1, i.e. Y satisfies the SDE 
dY t = (~f/2)Yf 1 + dB t , where B denotes a one-dimensional Brownian motion, the function 4>{x) = 
x 1_7 is harmonic for the process Y. A similar argument to the previous one permits to derive the 
attainability or non-attainability of the origin for the process Y, depending on the sign of 1 — 7. 

Roughly speaking, the strategy used in Delarue [I] to establish Theorem 13.41 aims to reduce the 
original limit absorption problem to the attainability problem for a suitable process. For this 
reason, we translate the absorption problem in Fq to the neighbourhood of zero. This amounts to 
investigate the behaviour of E(f e ), for £ e]0, 1], with X = (1, 1)* and % = inf{t > 0, X\ + Xj < £}. 
As guessed by the reader, our final result then writes in terms of a(0) and not of a(l, 1) as in Theorem 

m 

A crucial point is then to define the analogue of |-B| for the reflected diffusion X. In order to take 
into account the influence of the diffusion matrix a(0), we focus on the norm of X with respect to 
the scalar product induced by the inverse of the matrix o(0) (which is non-degenerate under the 
assumptions of Theorem 13. 4|) : 

Vt>0, Q t = {X u ar\G)X t ) 

= (1- s 2 r>rW + p 2 ~W -2 flP rV*M]- ' 

The process Q is devised to mimic the role played by \B\ 2 in the non-reflected Brownian case. Its 
differential form writes: 

Proposition 6.1 There exist a constant Cfe^J depending only on known parameters in Assumption 
(A), as well as a measurable function I j^"j| : M? — > M., bounded by Cfcj[ such that for N > 1 and 
t e [0,Cjv], C^v = mf{t > 0, \X t \ < N- 1 }: ' 

Vi G [0, Cat], dQ] /2 = T^X t )dt + ^Q^dt 

- WdH} + p^dHf] - Qt 12 [n\dKl + n\dKl\ (6.2) 

+ Q; 1/2 (a(X t )a- 1 (0)X t ,dB t ), 

with: 

k\ = [pf - sp^p^Xf] , K 2 t = [P2 2 ~ sp^P^X}} . (6.3) 

The proof of Proposition 16.11 is put aside for the moment. 

At this stage of the sketch, the most simple situation to focus on seems to be the case "s = 0" . In 
this framework, the dH term in ()6.2|) reduces to zero so that Q behaves like a standard Ito process 
in the neighbourhood of zero: 
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-1/2 

unc oucuai piuuutu 

: the process Q t 1 (a~ 

Brownian motion 



(1) For X close to 0, the scalar product Q t (a(X t )a l {Q)X t ,dB t ) driving the martingale part 



of Q 1 / 2 looks like the process Q t 1 (a 1 (0)Xt,dBt), which is, thanks to Levy's theorem, a 



(2) For X close to 0, the term I fej (-Xt) is negligible in front of Q i 



-1/2 



1 /2 

Thus, from (1) and (2), the differential form of the process Q t satisfies, for s = 0, the Bessel 
equation of index 1, at least for X in the neighbourhood of zero. 

For the sake of simplicity, we assume from now on that a(x) reduces for all x 6 [0, l] 2 to o(0), and 
in the same way, that I ju^ x) vanishes for x G [0, l] 2 . Up to the differential term dK, Q 1 / 2 then 
appears as a Bessel process. For this reason, we expect the function In to be a kind of Lyapunov 

1 /2 

function (in a sense to be precised) for Q t in the case s = 0. Ito's formula applied to (j6.2|) yields 
for Q different from zero: 

dln(Ql /2 ) = dB t - Qt 1 [K]dKl + K 2 dK 2 } . (6.4) 

It now remains to handle the dK term in (|6.4|) . Since s vanishes in this first analysis, the processes 
k 1 and n 2 driven by the differential terms dK 1 and dK 2 are always bounded from above and from 
below by positive constants. Recall moreover that Q^dK\ = Q^ 1 l^ x i =1 ydKl (the same holds for 
dK 2 ). We deduce that there exists a constant C > such that: 

dB t - C\dK\ + dK 2 } < d\n(Q] /2 ) < dB t - C~ l [dKl + dK 2 }. (6.5) 

The strategy then consists in comparing ()6.5|) to a bounded real functional of X. The most simple 
one is \X\ 2 . For b = 0, it is plain to derive from Ito's formula (note again that (Xt,dHt) = 
X^dHl + X 2 dH 2 vanishes since dH} = =0 ydHl and dH 2 = l^ X 2 =0 ydH 2 ): 



d\X\ 2 t = trace(a(0))dt - 2{X\dK\ + XfdKf) + 2(X t ,a(0)dB t ) 
= trace(a(0))cft - 2{dK} + dK 2 t ) + 2(X t ,a(0)dB t ). 

Add (|6.5jl and (|6.6|) (up to multiplicative constants) and deduce: 

d(-ln{Q] /2 ) + -^\X t \ 2 ) < -trace(a(0))dt - dB t + C(X t ,a(0)dB t ) 
d(-\n{Q] /2 ) + ^—\X t \ 2 ) > ( —tT&ce{a(<d))dt-dB t + C- l {X t ,a{<d)dB t ) 



(6.6) 



(6.7) 



1 /9 

Recall that Xq = (1,1)*. Take the expectation between and r^v = inf{t > 0, Q+ = l/A^}, 
for N > 1. Since the expectation of the martingale parts vanishes and since the matrix a(0) is 
non-degenerate, it comes for a new positive constant C'\ 

(C'y^HN) - 1) < E(rjv) < C"(ln(A0 + 1). 

Of course, E(tat) does match exactly E(T 1 /^r). However, it is well seen from the definition of Q 1 ^ 2 
that there exists a constant c > 0, such that t c -i^ < Ti/jy < t c at. This proves the third point in 
Theorem 13.41 □ 
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6.2 Non-Zero Cases 



As easily guessed by the reader, the cases "s < 0" and "s > 0" are more difficult. We derive first 
several straightforward consequences from Proposition 16,11 

(3) If s > 0, the dH term is always nonpositive. Up to a slight modification of the dK term (that 
does not play any role in the neighbourhood of the origin), the differential form of Q 1//2 writes 
as the differential form of Q 1 / 2 in the zero case plus a nonincreasing process. This explains 
why the process Q 1 / 2 reaches the neighbourhood of the origin in the case "s > 0" faster than 
in the case "s = 0" . We thus expect the case "s > 0" to be sub-logarithmic (again, in a sense 
to be precised). 

(4) On the opposite, if s < 0, the dH term is always nonnegative. With a similar argument to 
the previous one, we expect the case "s < 0" to be super-logarithmic. 

The strategy now consists in correcting Q 1 ! 2 to get rid of the dH term and to reduce the analysis to 
the one stated for "s = 0" . It is then rather natural to focus on the auxiliary process Vi > 0, At = 
Q] /2 + Z u with: 

vt > o, zt= -^=((p^, P ^)\x t ) = —L=[p^xl + P ^x 2 ]. 



s 



For s close to 0, the processes A and Q 1 ! 2 are equivalent in the following sense: there exists a 

constant C > 0, such that C~ l Q\^ 2 < At < CQ\^ 2 . The value of the constant C can be precised 

l li 

for s in the neighbourhood of zero. Indeed, the process Z can be expressed as: Z t = s^t{s)Q t > 
with 1 — e < 7t(s) < y/2 + e, for s in the neighbourhood of zero. 

Actually, we prove in Delarue |4_ that the equivalence property still holds true for all s g] — 1, 1[. As 
a first consequence, we derive that the processes Q 1 ^ 2 and A admit the same asymptotic behaviour. 
We thus focus on the second one and in particular on the differential form of A (the proof derives 
from Proposition 16. II and is thus left to the reader): 

Proposition 6.2 There exist a constant CfeT3 depending only on known parameters in (A), and 
a function TffiJ^ bounded by CftyJ^ such that: 

Vt e [0, Cat], dA t = IfcgX t )dt + ^Q; 1/2 dt 

- [R]dKl + R 2 dK 2 ] (6.8) 

(Qt 1/2 <j(X t )a-\G)Xt + -jJL^aiXtXpi^p^f^Bt), 

1 □ — 



with: 



Vt > 0, k\ = Qt 1/2 'n] + -jJ=pl l , 

- q; 1 ' 2 . 2 + ^L^p, 1 . 



Vi-s 

From now on, we assume that a and a do not depend on x, as done in the previous subsection 
to investigate the case "s = 0". In the same way, we assume for the sake of simplicity that r j|j^| 
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reduces to zero. Equation (|6.8|) becomes: 

Vi G [0,Cn], dA t = 1 + S ^ t{s) A^dt 



- [K\dKl + R 2 t dKj] (6.9) 



a, 



We try to mimic the arguments given in Subsection 16.11 For s close to zero, we deduce from 
Proposition 16. II that processes k 1 and k 2 are positive. Hence, from Proposition 16.21 the same holds 
for R 1 and R 2 . Since the general case s G] — 1, 1[ is really more difficult to handle, we restrict our 
analysis to s in the neighbourhood of zero. 

It now remains to investigate the martingale part in (|6,8|) . Its bracket writes: 

= [1 + 2Q; 1/2 Z t + ^((pj^p-y, aiO)^ 1 ^, 1 )*)] 
= [l + 2s 7t (s) + 2s 2 (l-s)- 1 ] 

Apply now the function x@ to A, for a given (3^0 (recall that this function is harmonic for the 
Bessel process of index 1 — (3). Ito's formula yields: 

dA p t = ^A^~ 2 Q t ((3, s)dt - 13 A^ 1 \R\dK} + R 2 t dK 2 ] + dM u 

with Q t ((3, s) = - 1)(1 + 2s 7t (s) + 2s 2 (l - s)" 1 ) + 1 + s lt (s)] 

where M denotes a martingale. 

We are now in position to complete the proof: 

1. If s > (and small), then the quantity @t{(3, s) tends to — s 7 t(s) — 2s 2 (l — s) _1 when (3 tends 
to zero. In particular, there exists f3 > (/3 close to zero) such that @t((3, s) < 0. For this 0, 
dA^ < -j3Al~\R\dK\ + R 2 dK 2 } + dM t . We thus recover the right hand side in (jB3)) with 
respect to the function x@ and to the process A. We then derive the second inequality in 
(|6.7|) . with — ln(x) replaced by — x@, and Q x l 2 by A. 

This proves the first point in Theorem 13.41 

2. If s < (and small), we can follow the previous argument and thus prove that there exists (3 < 
{13 close to zero) such that @ t (f3, s) > 0. For this 0, dA P t < -l3A P t ~ l {R\dK} + R 2 t dK 2 } + dM t 



Again, we recover the right hand side in (|6.5fl with respect to the function x@ and to the 
process A, but with — C _1 > 0, i.e. C < 0. We then derive the second inequality in 1)6.7(1 . 
with — ln(x) replaced by — x@, Q 1 ^ 2 by A and C by —C. 

This proves the lower bound in the second point of Theorem 13.41 (with respect to a possibly 
different (3-). 

To obtain the upper bound, note for (3 — > — oo that Ot(/3, s) ~ (3{1 + 2s 7t (s) + 2s 2 (l — s)^ 1 ) 
and deduce (at least for s close to zero) that there exists f3 < such that @t((3, s) < 0. For 
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this (3, dA t > —f3At [R\dKj + R 2 dK 2 ] + dM t . Again, we follow the argument given in 
the case "s = 0" and then derive the upper bound in the second point of Theorem 13.41 (the 
resulting /3 + may differ from the one given in the statement). □ 



6.3 Proof of Proposition I6TT1 

Write first the SDE satisfied by the process Q: 

Lemma 6.3 Keep the notations introduced in the statement of Proposition Wl\ Then, there exists 
a constant GtffTfo depending only on known parameters in (A), as well as a measurable function 
Ijg^j, bounded by Offi^ such that: 



Vt > 0, dQ t = Ql /2 IftjiX t )dt + 2dt 

^^Q\ 12 [Pi l dH l t + p^dH 2 ] - 2 [nldK] + K 2 dK 2 ] (6.10) 

+ 2{a(X t )d- 1 (0)X t ,dB t ). 

Assume for the moment that Lemma 16.31 holds and complete the proof of Proposition 16.11 Apply 
Ito's formula to Q 1 ! 2 and derive from the differential form given in Ijfi.lOJ) (note that Q does not 
vanish for t G [0, Cn])'- 



Vie [0,Cat], dQ] /2 = -JfegXJdt + Q, 12 



dt 

[p^ l dH l t + p^dH 2 ] - Q; 1/2 [K\dK l t + K 2 dK 2 } 



Vi^V 2 

+ Q; 1/2 (a(X t )d- 1 (0)X t ,dB t ) (6.11) 
-^Qf ,2 (X u d-\Q)d(X t )d- l ($)X t )dt 

= [A(l,X t )+Q; 1/2 ]dt + dA t (2)+dA t (3) - ^Q; 3/2 A(A,X t )dt. 

Focus on A(4, •) in (|6~TT|) : 

Vt e [0,(n], A(A,X t ) = Q t + (X^arHoMXt) - a(0)]a- 1 (0)X, 
= Q t + A(5,X t ). 



(6.12) 



Thanks to Assumption (A) (boundedness, ellipticity and Lipschitz continuity of a), the function 
A(5, •) satisfies: 

Vt € [0,0v], \A(5,X t )\ < CQ 3 / 2 . (6.13) 
where the constant C depends only on known parameters. Plug now (|6.12l) into (|6.11|) : 

Vt G [0,0v], dQ 1 / 2 = [A(l,X t ) - lQ; 3/2 A(5,X t ) + ^Q; 1/2 ]dt + dA t (2) + dA t (3). (6.14) 



Set finally IfojjXt) = A(l,X t ) - ^Q t 3/2 A(5,X t ). Thanks to (l6~TT|) . (l6~T3l and to (l6~T4l . this 
completes the proof. □ 
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Proof of Lemma 16.31 Thanks to Ito's formula, d(X,a 1 (0)X) writes: 

Vt > 0, d{X t ,a- l {V)X t ) = 2(X t ,d~ 1 (0)b(X t ))dt + tvajce[a- 1 (0)a(X t )a- 1 (Q)]dt 

+ 2{X t ,d~ 1 (0)(dH t -dK t )) + 2{d-\0)X t , a(X t )dB t ) 
= 2dt (6.15) 
+ 2(X t ,d- 1 (0)b(X t ))dt + trace [CT~ 1 (0)(a(Xt) - a(0))a _1 (0)]rft 
+ 2<X t ,a- 1 (0)(dff t - dK t )) + 2(a- 1 (0)-X t ,?(X t )dB t ) 

Set J^gX t ) = Q t - 1/2 [2(X t ,a- 1 (0)6(X t ))+trace[^- 1 (0)(a(X t )-a(0))a- 1 (0)]]. 
Focus now on the <iff and dif terms in the above r.h.s. Due to (|3.3|) . note first that X^dH^ = 
and X 2 dH 2 = and that p^XfdH} = (1 - s^^Q]' 2 dH} and p\ x X\dH 2 = (1 - s 2 ) l l 2 Q l t /2 dH 2 . 
Derive that: 

(Xt.g-Hojdfft) = ^(a-HO))^* 1 +X|(a- 1 (0)) 22 dif t 2 
+ (a- 1 (0)) lj2 [^F f 1 +X t 1 ^ t 2 ] 

= -r^r V + *^f] (fU6) 



g{ /2 [ P r 1 dff t 1 + A'2 1 dff*]- 



s 



Note again from JH3J) that X/diT/ = dK\ and X t 2 <iFQ 2 = dK" 2 . Thus: 

{J^tT^O)^*) = X^a-^O))^^ 1 + Xf(a- 1 (0)) 22 dK 2 

+ (a-^O))^^ 1 ^^^ 2 ^ 1 ] (61?) 

= [/°r 2 - ^rVj 1 ^ ] ^ + 3-^2 [P2 2 - »PF W*t] dKl 

Plug now 16|) and (|6.17() into ()6.15f> and complete the proof of (|6.10f) . □ 

Note from (IfHTI) . iKTTSll . (lfi~TH) and (ffTTBl that Ijr^i reduces to zero if b vanishes and a is constant. 
This explains why we can assume in Subsections ^ and !6.2l d(x) = a(0), 6 = and f]|^= H6~2] = 
0. 



7 Conclusion 

Probabilistic models depending on an extra random environment have been widely developed for 
twenty years. From a practical point of view, the environment expresses external constraints that 
satisfy some statistical properties. For example, the behaviour of an investor that borrows a cer- 
tain amount of money may depend on the values of various assets as well as some quasi-periodic 
parameters. 

From a more theoretical point of view, our work shows that earlier results (see e.g. Louchard |17| . 
Louchard et al. |18j and |19j and Guillotin and Schott 10 ) for the diffusive behaviour of the banker 
algorithm extend to a more general framework. This is a new small step in the robustness of the 
underlying model. The theory of stochastic homogenization is also relevant for other distributed 
algorithms such as the problem of colliding stacks. 
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The reader may object that our own paper exclusively focuses on the trend-free case (see Assump- 
tion (A.4)), whereas so-called contracting and expanding regimes for distributed algorithms have 
been investigated in the previous references. We reasonably guess that the strategy presented here 
applies to the expanding case, after a suitable substraction of the trend, as done for example in 
Louchard and Schott 18 . On the opposite, the contracting framework would require further devel- 
opments, both based on homogenization techniques and on large deviations arguments as in Maier 
|20| . We leave this for further research. 

In a different perspective, two questions follow from the study of the limit reflected diffusion. On 
the one hand, our global result writes as a double asymptotic property: we let first m tend to 
+00 in m~ 2 T( m ) and then A to 2 in E(T\) (at least in the two dimensional setting). A natural 
question would consist in investigating m _2 E(T^ ) (with obvious notations) as m tends to +00 
and for (m _1 A m ) m >i converging towards 2. On the other hand, the extension of Theorem 13.41 to 
the d-dimensional case appears as a challenging problem. Nevertheless, numerical methods offer 
a possible alternative to obtain an empirical description of the asymptotic behaviour of the limit 
diffusion or even to compute with a Monte-Carlo procedure an estimate of underlying quantities of 
interest. 
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